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Abstract 

We present analytical solutions to the steady state injection-condensation-coagulation equation of 
aerosols in the atmosphere. These solutions are appropriate under different limits but more general 
than previously derived analytical solutions. For example, we provide an analytic solution to the 
coagulation limit plus a condensation correction. Our solutions are then compared with numerical 
results. We show that the solutions can be used to estimate the sensitivity of the cloud condensation 
nuclei number density to the nucleation rate of small condensation nuclei and to changes in the 
formation rate of sulfuric acid. 
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1 Introduction 


Aerosols in Earth’s atmosphere are important for the climate systems as they are required to form clouds. 
As a consequence, different aerosol characteristics translate into different cloud properties, such as dif¬ 
ferent radiative forcing, different life time and precipitation. Since a large part of the characteristics of 
aerosols can be described by their size distribution, solving for the distribution is of particular interest. 

For a spatially homogeneous aerosol distribution, the size distribution can be described with a density 
function n{v, t) where v is the aerosol size, such that n{v, t)dv is the total number of aerosols per 
unit volume, with volumes between tii and V 2 - 

The main equation describing the temporal evolution of n{v, t) is the coagulation-condensation equation 
(e.g., Peterson et al, 1978; Seinfeld & Pandis, 2006): 

dll d 1 

^ (■y, f) n] = -J f3 {v - v,v)n{v - v,t)n{v,t) dv 

(3 {v, v) n {v, t) n {v, t) dv — R{v, t) + S{v, t), (1) 

where I{v, t) = dv/dt is the growth rate of a particle of volume v, /3(ni, ^ 2 ) is the coagulation coefficient 
for particles with volumes vi and V 2 , R{v, t) is the rate with which particles are removed from the system, 
and S{v, t) is the nucleation rate of fresh particles. It is of course similar to the Smoluchowski (1916) 
equation, with the addition of the condensation term. 

Klett (1975) provided an analytical solution to the above equation, but without the condensation term, 
for a few specific cases by means of a Laplace transform. For example, he solved the aerosol distribution 
for the case of a mass independent coagulation coefficient. 

Later, Ramabhadran et al. (1976) provided an analytic solution to the more general problem which 
includes condensation. Like Klett, they have shown that power laws can provide approximate solutions 
in different regimes. Peterson et al. (1978) elaborated on the above and provided more general solutions. 
For example, they solved the time dependent problem. However, they confined themselves to cases in 
which the condensation is either constant or linearly dependent on the aerosol volume. 

It should be noted that an equation similar to the aerosol condensation/coagulation equation appears in 
other contexts as well. First, the original coagulation equation (with no condensation) was formulated by 
Smoluchowski (1916) to describe a colloidal fluid. A second example is that of cloud raindrops which 
behave similarly to the smaller sized aerosols with the main exception that their coagulation equation 
includes raindrop fragmentation (e.g., see review by Beheng, 2010), while outside the terrestrial settings, 
Birnstiel et al. (2011) solve a coagulation/fragmentation equation describing dust grains in the interstellar 
medium. 

We begin in §2 by writing the coagulation-condensation equation in a dimensionless form and the ap¬ 
proximations we assume in the present analysis. This will define the problem that we solve. In §3, we 
first solve the problem while discarding the condensation term, and by doing so we arrive at the same 
solution found in Klett (1975) by different means. This solution is expanded by treating the condensation 
term as a perturbation. The new analytic solution is appropriate for large aerosols where the dominant 
process is coagulation. In §4, we solve the opposite limit, first, when the coagulation term is altogether 
negligible, and then, when its effects are approximated. This solution is appropriate for small aerosols 
where the dominant process is condensation. In §5 we compare our results to a numerical calculation 
(the details of which can be found in the appendix). We end in §6 by discussing the significance of our 
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theoretical results. 


2 The problem 

Although the general problem of aerosol growth, eq. 1, is time dependent. We will concentrate in the 
present work on the time-independent case; we will assume a steady-state solution n{v). 

As in Peterson etal. (1978), we will assume a constant coagulation coefficient /3 = /3o and / {v, t) = av'^ 
where 0 < 7 < 1 , and we will emphasize the case 7 = 1/3 (which is appropriate for the “continuum 
regime”) and 7 = 2/3 (which is appropriate for the “kinetic regime”). The continuum and kinetic 
regimes are defined by the Knudsen number Kn = X/R, where A is the mean free path, and R is the 
radius of the paricle. The continuum regime corresponds to Kn <C 1, and the kinetic regime corresponds 
to Kn :§> 1. The transition between these two regimes actually takes place for Kn « 0.1 (at which point 
the interaction rate is half the kinetic regime), occurs when R ^ X ^ 0.5/um (for the earth’s atmosphere 
at sea level, e.g., see Seinfeld & Pandis 2006) 

We also assume that the timescale for removal of particles from the system is much larger than the time 
necessary to arrive at the steady-state solution. Therefore we will discard the removal term altogether. 

Finally, we assume a delta function source S {v,t) = S06 {v — vi), again following Peterson et al. 
(1978). As we shall see below, we will also assume that is much smaller than some characteristic 
volume that will be defined lafer. 

We fherefore obfain fhe following inlegro-differenlial equafion for fhe sfeady-sfafe disfribufion: 

d I r /■“ 

— [av'^n {v)] = - / Pan (v) n (v — v) dv — / fio'iT'(v) n (v) dv + So5 {v — vi). (2) 

dv 2 Jq Jq 


2.1 General Kernel 

Allhough we solve in the present work the coagulation equation while assuming a constant coagula¬ 
tion coefficient, we note that the any solution under this assumption can be immediately generalized to 
coagulations having kernels of the form /3 {v, v) = 

Here, the coagulation equation has the form 


A 

dv 


1 rv /•oo 

[crv'^n (v)] = - / f3iv°‘(v — v)°'n{v)n {v — v) dv — / /3iv‘^v‘^n{v)n (v) dv + SqS {v 

2 Jo Jo 


To reduce it to the fixed kernel form, we define 

f{v) = v°‘n{v), 


t>i). 

(3) 

(4) 


giving rise lo fhe following equafion for /: 
d ^ 

— [cru^-"/(u)] = - y /3if{v)f{v-v)dv-J I3if{v)f {v)dv + Sod {v - vi) . (5) 

This equation is fhe same as eq. 2, bul wilh a differenl value of 7 , fhaf is 7 —)> 7 — a. 
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2.2 Dimensionless Form 


Before we solve eq. 2, we recast it in a dimensionless form. We begin by defining xo as the total number 
of particles in the system, that is, xo = n (v) dv. The differential equation which describes its 
temporal dependence is (Klett, 1975): 


^^0 _ ^ « 2 I c 


( 6 ) 


which gives a steady state solution of xo = y^25o//3o- 
Next, we define fhe characferisfic volume as 


^2 \ 1 /( 2 - 27 ) 


V2 = 


and characferisfic particle density fo be 


no 


(^) 

y/25o//3o 


V2 


(7) 


( 8 ) 


We also defined a dimensionless lime as t = t/y/2 /SqI3q. 

Using Ihese characlerislics quantifies, we can define fhe dimensionless variables as x = v/v 2 , xi = 
v\lv 2 , and y = n/no, such fhaf fhe time dependenf version of eq. 2 becomes 


dr 


d 

+ \/2— [x^y (x)] = / y {x)y {x — x) dx — 2y (x) + <5 (x — xi). 
dx Jq 

Nofe fhe facl fhaf in sfeady sfafe: 


y (x) dx = 1. 


So fhe sfeady slate equation can also be written in fhe form 

[x'^y (x)] = / y (x) y {x — x) dx — 2y (x) / y (x) dx + 6 {x — xi) 
dx Jo Jo 


(9) 


( 10 ) 


( 11 ) 


which is closer lo fhe form of equation 2. We will assume Ihroughouf Ihis paper fhaf xi <C ^2 , which 
Iranslafes lo xi <C 1 in fhe dimensionless quantifies. 

The amounf of condensable matter is affecled by fhe processes of nuclealion and condensalion. Mosl of 
Ihis matter is sulphuric acid (S A) and water which condenses wilh if (Ihough over fhe oceans melhanesul- 
fonic acid, MSA, could also be importanf). The equation which describes fhe change of fhe condensable 
malerial is: 


dMsA 
dt 


— ^m. 


/*o 

J Vl 


' Q 

(TV ’ —dv — bovi- 

dv 


( 12 ) 


where em is fhe rale of creation of sulphuric acid (or MSA) by an exlernal source, logelher wilh fhe wafer 
which Icondenses wilh if. Therefore, in sleady-slale we have: 


dv 

em= I av'^—dv - SoVi. 


(13) 


ivi 


Laler on we will discard fhe SoVi ferm as if is smaller fhan fhe firs! term (under fhe assumpfion vi <C X 2 ). 
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3 The coalescence limit 


3.1 Rough approximation 

A full solution to the coalescence equation will be presented in §3.3. However, before we do so, it is 
worthwhile to derive the general behavior from very simple arguments. This will not provide a solution 
with accurate normalization constants, but it encapsulates the underlying physics and therefore it provides 
the correct power law. 

We start with eq. 2, and define {v) to be the rate of change in the total volume of particles of size v or 
larger, that is, 

roo Q 

Cm (v) = J —n{v)vdv. (14) 

We shall try to calculate this quantity crudely, and then require it to be independent of v, as we expect it 
to be in steady state. We shall first assume no condensation, and later generalize to the case also having 
condensation as a perturbation. 



Figure 1: is estimated by adding the contribution from the shaded area in the (wi,« 2 ) plane—the 

volumes of the coagulating aerosols. Unshaded regions don’t contribute to em either because they don’t 
generate more volume above v, or because there are no particles in those regions (with ui or U2 < v. In 
each region, the contribution to the integral behaves differently and therefore its typical value is different, 
as is summarized in table 1. The light shaded area is a mirror of the darker shaded one. Instead of 
counting it separately, one can simply double the contribution from the darker region. 

The contribution to em from coagulation comes from a double integral over pairs of aerosols {ui,U 2 )- 
Without loss of generality, we assume ui > U2, and we only need to consider cases in which U2 < v but 
ui + U2 > V, because coagulation between two aerosols larger than v contributes nothing to em- We can 
then distinguish between 6 different regions in the remaining part of the (ui, M 2 ) plane, as is depicted in 
fig. 1 . The approximafe contribution from each of these regions to em is summarized in table 1 . 

For example, we shall consider the case where the volumes of the coagulating of both coagulating parti¬ 
cles is smaller but of the same order as v. As a consequence, the volume “produced” from the coagulation 
is V, the area in the {ui,U2) plane is of order Since both coagulating particles are of the same order 
as V, their total number will be of order n{v)‘^, i.e., we assume it to be constant in the integral, which is 
why the total contribution from this region comes out to be ~ /3o x x n{v)‘^ x v = / 3 ov^n{v)‘^. 
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Number 

Definifion 

Area 

Number of 

Parficles 

Volume 

confribufed 

Tofal confribufion 

to €m 

1 

^ 1,2 ~ v 

v'^ 

n{vY 

V 

(3QV‘^n{vY 

2 

Ul ^ V, U2 ^ V 


n{v)‘^ 

V 

l3QV^n{vY 

3 

Ul ^ V, U2 ^ V 

00 

n{v)n{ui) 

V 

Pov‘^n{v) f^n(ui)dui 

4 

Ul ^ V — Vl, U2 ~ 


n{v)n{vi) 

V 

/3ovfvn(v)n(vi) 

5 

Ul '^V,U2% f 1 

VlV 

n{v)n{vi) 

Vl 

/3ovfvn(v)n(vi) 

6 

ui > n, M2 ~ '^’1 

00 

n{vi)n{ui) 

Vl 

/3ov^n(vi) f^n(ui)dui 


Table 1: Regions in the {ui,U 2 ) plane contributing to e^. For each of the regions in fig. 1, we crudely 
approximate (up to numerical factors) the area of the region in the {ui,U 2 ) plane, the number of co¬ 
agulating particles n{ui)n{u 2 ), and the total volume contribution to particles of size v or more. Using 
these quantities, we estimate the total contribution to as /3oX Area x number of particles x volume 
contributed. We deal with infinities by later assuming a power law behavior for n{v) and calculate the 
integral. 

We now assume a power-law solution, n{v) = Av~'p. It is then apparent that one must have p > 1, 
otherwise the total number of aerosols f n(v}dv will diverge. We also assume p < 2, otherwise the 
integral f n(v)vdv converges to a finite volume and not one which can increase linearly with time, as we 
expect in the steady state. 

Under the aforementioned assumptions to total contribution to Cm from the regions described in table 1) 
comes out to be 

em(v) ^ /3ov^A^v~^P -h /3ovfvA^vf^v~P ^ /3 oA^v^~^p, (15) 

where we have also assumed vi v for the second approximation. 

If we further assume that em(v) is independent of v, we find fhaf p = 3/2 and A ss s/^mJWtt- As we 
shall see in fhe more rigorous freafmenf described in §3.2, fhis resulf is correcf up fo a numerical facfor. 

The lasf resulf can be furfher improved by adding fhe leading correcfion obfained when including con- 
densafion. To do so, we guess a solufion in fhe form niv) = Av~'p + Cv~^. This leads fo fwo addifional 
ferms in fhe expression for Cm, which are 

poo 

Cmiv) ^ (3oA^v^~^^ + 2/3oACv^~^~‘^ + / Au~^av?du + Av^~'^av'^, (16) 

J V 

where fhe second term is a second-order coagulafion ferm, fhe fhird ferm is fhe firsl-order condensation 
term describing fhe growfh of particles larger fhan v, while fhe lasf ferm describes fhe flux of particles 
which fhrough condensafion become larger fhan v. Nofe fhaf fhe lasf fwo ferms are similar in size if fhe 
fhird ferm converges which we require anyway. This requiremenf is equivalenf fo 7 < 1/2. Wifhouf 
fhis requiremenf, an upper cutoff for fhe disfribufion is necessary fo limif fhe condensafion onfo fhe large 
volume fail. 

The “zerofh” order equafion obfained from requiring fhat musf be independenf of v in sfeady sfafe 
leads fo fhe same values of A and p found above. The “firsf” order equation obfained leads fo 

0 « (17) 

Since fhis equation holds for all v, we find fhaf q = 2 — 7 and C ^ —ajfdQ. As we shall see in §3.3, a 
more rigorous freafmenf yields fhe same correcfion up fo a 7 -dependenf numerical facfor. 
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3.2 Full Solution with no condensation correction 


Klett (1975) solved the fixed coalescence cross-section case without while neglecting the condensation 
equation using Laplace transform, and found that 


n (v) 



(18) 


for large values of v. Here we show a shorter and more intuitive way to reach the same solution. In 
the next subsection, we will add the correction term obtained when adding the condensation at large v’s 
(when the condensation term is necessarily small). 

We define e (x) to be the dimensionless rate of change of the total dimensionless volume of particles of 
dimensionless size x or more (with “dimensionless” hereafter omitted), that is: 


e(x) = 



—y(x)xdx. 

OT 


(19) 


Since we neglect condensation, this increase in the total volume is due to two “types” of particle coa¬ 
lescence. First, two particles with a volume smaller than x can combine into one particle with a volume 
larger than x. Second, a particle with a volume larger than x can combine with a particle with a volume 
smaller than x. Together, we find that the total volume change due to coalescence is 

f'X f*X /*oo 

^{x)= / y (6) y (6) (6 + 6) + 2 / / y (6)y (20) 

7^1=0 7^2=3:—?l 7^1=0 7^2=3: 

which simplifies to 

e{x) = 2 / y (?i)y (21) 

76=0 

Next, we assume a power law solution, i.e., y = Bx~^, and that 1 < p < 2. The assumption p > 1 is 
necessary to ensure a finite total number of particles, while the assumption p < 2 is necessary to ensure 
an infinite total volume. The latter is necessary because the total volume is a monotonically increasing 
function of time. Together with eq. 21, we find affer some algebra that 

e{x) = 2 B‘^x^~‘^^ - / {1 — dr], (22) 

P - 1 Jri=0 


In steady state, e should be independent of x. Therefore p = 3/2, and 

e = 2B^—^[^ {1 - drj = AttB^ (23) 

oj i — 7^=0 

Thus, the steady-state solution to eq. 9 is 



If we define e™, to be the rate of change of volume in the original equation, then it can easily be seen that 

Em = 5'of2e- (25) 


7 






Therefore, the solution to eq. 2 with no condensation (for large aerosol volume) is 


which is the same as the solution found by Klett (1975). 
Incidentally, we can use equation 25 and the definition of V 2 to get 

,_ / p \ (1-7) 


This result will be used later in order to eliminate a from our results, by introducing the dimensionless e 
whose value depends on 7 . This is useful because under most physical scenarios, cr is unknown, but 
is determined from various chemical or physical processes. 


3.3 Full Solution with condensation correction 


The next step is to generalize the calculation described in §3.2 by adding the first-order correction term 
associated with condensation. With the latter term, the equation describing the increase of the total mass 
of particles larger than x now includes two additional terms 

e{x) = 2 / y i^i) y {^ 2 ) ^idCidC 2 + V2Cy{0dC + ^x'^y{x)x. (28) 

7 ^ 1=0 7 ^ 2 = 3 ;—Jx 


The second term in the equation describes volume change due to condensation on particles of volume x 
or larger. The third term describes particles of volume slightly less than x, which in a unit time grow to 
become larger than x due to condensation. 

We will now look for a solution of the type y = Bx~p + Dx~^, where q > p. We now plug this solution 
into the equation for e, and neglect small powers of x. For 7 < 1/2, the equation for the highest power 
of X is unaffected by the new term Dx~'^, so we still have p = "ij2 and B = e/( 47 r) as before. 

The equation we get for the next highest power of x is (after some algebra): 


0 = 2Dx^-P-'^ 


1 1 

-7 7-7 

q — I p — 1 



{1 - r])^-Ur] + V2x^+'^-P 


1 + 


p- 1 -7 


(29) 


As before, this equation should hold for all values of x (in the limit x ^ 1). Equating the power law 
indices of x gives us: 

2-7 = 9- (30) 

This explains why 7 should be smaller than 1/2 for this solution to be valid. Next, we can equate the 
coefficients. Once we plug in p = 3/2 and q = 2 — 7 , we find 


2D 


1 1 
+ 


[1-7 1/2 J 


^ 1/2 _ y^^7 ^ d-q = —\f2 ( 1 + 


1 


1 / 2-7 


(31) 


This equafion can be solved for D. For example, in the important case where 7 = 1/3 we have 
/q 7y“^/^(l — p)3“^dp = yTrr(l/3) /T (5/6) Rs 4.20655 such that eq. 31 gives 79(7 = 1/3) = 
—0.3362. For other values of 7 , see fig. 2. 
















Figure 2: D vs. 7 for the case 1/2 > 7 > 0. 


To summarize, the solution of eq. 9 is of the form (for 0 < 7 < 1/2): 


y{x) 



3/2 +^^-( 2 - 7 ), 


For the specific case 7 = g, we obtain 


y{x) 



0.336x"^/^. 


Using the physical quantities, this specific solution becomes 


n (v) 



- 0 . 475 ^?; 
Po 


■ 5/3 


(32) 


(33) 


(34) 


In many cases, we are given Cm, i-e., the sulfuric acid formation rate, instead of a, its equilibrium number 
density. We can therefore use eq. 27 and find 


n (v) 



0.475^o“^/^5o“^/® 



(35) 


A numerical solution of the equation (described in the appendix) shows that e (7 = 1/3) Ri 3.296. This 
is perhaps the most important result in the present work since it describes the leading two terms in the 
distribution of aerosols in the atmosphere. The analytic form of the second term is described here for the 
first time. 

In principle, it is possible to add higher order corrections as a power series of x, where each term’s power 
decreases by 1/6. 


3.4 Full Solution with condensation correction for 7 > 1/2 

The solution described in §3.3 is valid only for 7 < 1/2, because it was assumed that q > p such that 
2 — 7 > 3/2 (using eq. 30). Therefore, the larger 7 case should be solved separately. As we shall see 
below, this has a major effect on the solution—the leading power of x changes. 
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First, because e{x) diverges, instead of working with e{x) we will work with e(x, X 3 ) which is the rate of 
change in the total mass between x and X 3 . In steady state, it clearly vanishes, even though the integral 
above X 3 diverges. 

Again we assume a power-law solution of the type y = Bx~^, and write the solution in the form 

e{x,X 3 ) = F{x) - F{x 3 ), (36) 

where the expression for F{x) is exactly the same as the expression we had for e(x) in eq. 28, where we 
drop the upper bound wherever it diverges. Using this form, we find after integration that 

F{x) = [ rj^-P (1 - rj)^~P drj + V2Bx^+^-P - V2B - ^(37) 

P - 1 7 , 7=0 1 + 7 - P 


Since e(x, X3) should vanish, we must require F{x) to be independent of x. However, since p / 3/2, 
the only way for F{x) to be independent of x is for the x-dependent terms to cancel each other out. This 
implies that the exponents are equal, such that, 2 — 7 = p. Last, we can find the pre-factor B by requiring 
the coefficient of x^“^^ (or x^^'^~P) to vanish as well. The result is 


y (x) = Bx^ 


/o ~ 7 )^ 1 


(38) 


For instance, the solution to eq. 9 for 7 = 2/3 and large values of x is approximately given by y = 
0.230x“‘^/^. In terms of v, the solution to eq. 2 for the above case is 

n{v) = (39) 

Po 


For other values of 7 see fig. 4. Unlike the previous case, of 7 < 1 / 2 , the condensation term’s contri¬ 
bution to em formally diverges. This implies that without additional physics introducing a large volume 
cutoff, such as dry or wet deposition, or a dependence of 7 on the aerosol volume, the solution is no 
longer physical. It also implies that we cannot eliminate a from the above results. 


3.5 Full Solution with 2"*^ condensation correction for 7 > 1/2 


As for the small 7 case described in §3.3, it is possible to derive a higher order correction to the lowest 
order solution described in §3.4 above. To do so, we can repeat §3.4 under the assumption that the 
solution is of the more general form y = Bx~p + Dx~p (where q > p), and calculate F{x) defined in 
eq. 36. 

We then require F{x) =0 and compare the leading terms in powers of x. The lowest order term was 
described above (see eq. 38). The next order gives. 


/■i 1 

0 = 2BDx^-P-'^ / - (l-rjV-Ur] + 

Jo 9-1 

/■i 1 

+2BDx^-P-^ / r,^-^ - {l-r]f-Pdr] + 

Jo P - 1 


+V2Dx^^'^-'^ - V2D 


1 


1 + 7-9 




(40) 


10 










The requirement on the exponents gives ‘i—p — q = l + 'y — q, however, it follows straightforwardly from 
p = 2 — 7 found in the first order solution. Consequently, we only obtain new information by requiring 
the sum of the coefficients to vanish. Since D cancels out, the constraint becomes 




1 

1 + 7-q 


V2 = 2B 


1 1 

-7 “I-7 

q — I p — 1 



(1 — r])^ dp. 


(41) 


This equation can be solved numerically for q, given the previously calculated values of p and B. For 
example, the solution for the case 7 = 2/3 is gsa 1.47. Solutions for other values of 7 are described in 
fig- 3. 



Figure 3: q vs. 7 for the case 1 > 7 > 1/2. 



Figure 4: B vs. 7 for the case 1 > 7 > 1/2. 

The value of D is found numerically, and it is found to have a rather strong dependence on xi. Note that 
D can be negative. The numerical simulations also seem to indicate that for xi —0, D{xi) —)• —B. 
This can be seen in fig. 5. 
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Figure 5: (color onloine) D vs. 7 and —B vs. 7 for the case 1 > 7 > 1/2. The dashed line depicts —B 
while the solid lines are from top to bottom: xi = 10,1,0.1,0.01,0.001. 


In terms of the equation 2, we have the following solution for the case 7 = 2/3: 


n {v) = 0.325^u + D{xi = vi/v 2 )^j^^ ^ 


a 


Sof3o) 


0.706 


-1.471 


(42) 


4 The condensation limit 

Until now we studied the limits where coagulation is much more important than condensation. We 
now concentrate on the opposite limit, where condensation is much more important. This describes, for 
example, the growth of small aerosols. We begin by describing the pure condensation limit, and then 
continue by crudely adding the coagulation as a correction. 


4.1 The condensation solution without coagulation corrections 


After discarding both coagulation terms in eq. 9, the steady state equation can be solved analytically. Eq. 
9 simply becomes 


\/2— [x'^y (x)] = 6{x — xi). 
ax 


(43) 


The general solution is y 
yields 


Cx where C can be found by the “boundary condition” at x 


y = Cx . 

V2 


xi, which 
(44) 


4.2 The condensation solution with a coagulation correction 

The pure condensation solution can be extended by adding the first order correction arising from coag¬ 
ulation. This introduces two terms appearing in eq. 9 which are estimated by using the “zeroth order” 
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solution. The first coagulation term is roughly y {x) y {x — x) dx ^ y{x)y{x)x ^ x^~‘^'^. The second 
coagulation term is roughly —2y{x) « x~'^ . Therefore, assuming 7 < 1, the first coagulation term will 
be much smaller than the second coagulation term for all a: <C 1 , and we shall not take it into account 
in the calculation that follows. The steady state dimensionless equation can still be solved analytically. 
The equation is 

V2-^[x'^y{x)] = -2y{x) + 5{x-x\). (45) 

Its general solution is given by 


y = Cx exp 


_V2— 

1-7 


Using the “boundary condition” at x = xi, and assuming xi <C 1, we obtain: 


(46) 


y{x) 


^x--exp 


_V2—x'-^ 

1-7 


(47) 


or in terms of eq. 2 

cr 1 — 7 

In the case 7 < 1/2, we can eliminate cr from this result, and give it instead in terms of em, 
This gives 



a 


n{v) « So/(do 



V exp 



(48) 

using eq. 27. 

(49) 


where e is 7 -dependent, and must be calculated numerically. 


Note that this solution can easily be generalized to the case where a particle-loss term — An (x) is added 
to eq. 2 . 


5 Comparison with the numerical solutions 

Although there is no exact analytic solution for the full coagulation-condensation equation, a full solution 
can be obtained numerically, as is described in the Appendix. 

The full solution can then be used to check the quality of the analytical approximations obtained above. 
Fig. 6 plots the full numerical solution and the analytical approximations obtained for the case 7 = 1/3. 
As is evident, the two analytic approximations for the x 3 > 1 and the x <C 1 limits are relatively accurate. 
In fact, at X ~ 0.5 both solutions are only 30% off the exact solution. 


6 Discussion 

One of the interesting implications of the solution to the coagulation/condensation equation is that it 
helps us understand how the aerosol distribution will change under different perturbations to the CN 
nucleation rate and to the amount of available condensable gas. 

It is known from numerical simulations that perturbing the nucleation rate So will have an appreciable 
effect on the density of small CNs, whose growth is governed by condensation. However, the effect is 
going to be small for large CCNs once coagulation becomes important (Svensmark et ai, 2013). 
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Figure 6: (color online) Comparison between the numerical solution (blue) and the analytic approxima¬ 
tions derived here, for the case 7 = 1/3. The zeroth order solution describing coagulation and its first 
order correction from condensation are depicted with the red and cyan lines respectively, which are valid 
for large x’s. The zeroth order condensation solution with its first order coagulation correction valid for 
small x’s is described by the green line. 


On the other hand, a perturbation to the amount of available condensable material em, which is usually 
sulfuric acid, will affect both the number density of small CNs and of the larger CCNs. 

With the result of eq. 35 and eq. 48, this behavior can be understood and quantified analyfically. 

If is probably possible fo obfain beffer fheorefical solufions in fhe coagulafion limif by faking fhe nexf 
ferm in fhe power series. However, this will probably not introduce any new physical consequences to 
the solution—unlike the first order correction to the solution (see eq. 35) which added a dependence on 
the nucleation rate. 

We also solved the case 7 > 1/2, which corresponds to high condensation rates for large particles (see 
eq. 38). The zeroth order solution was again improved by finding fhe firsf correction ferm. Since fhe rafe 
of change in fhe fofal volume diverges, fhe solufion has no physical meaning wifhouf infroducing, for 
example, a large size cufoff. 

As menfioned in §2.1, anofher inferesfing aspecf is fhaf fhe solufion can be sfraighfforwardly generalized 
fo describe kernels of fhe fype /3 (n, v) = fhrough fhe definition 7' = 7 — a. This allows us fo 

describe more realistic aerosol growfh in fhe real afmosphere. 

This implies fhaf in fhe condensafion limif fhe generalized solufion is 


y{x) 


v/2 


— w— 2ck 

X ’ exp 


-v/2-^-X 

1 — 7 — a 


1—7—a 


(50) 
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while for the coagulation limit we find 


y{x) = Bx ^ ^ + Dx ^ (51) 

where B,D,p and q are defined with 7' instead of 7. 

For example, In the limit were coagulation is through hydrodynamic capture of particles in the Stokes 
regime (Klett, 1975), we can expect 0 = 2/3. 
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Appendix - Numerical Solution 

Here we describe the numerical solution of eq. 2. Since the equation is the time-independent limit of the 
full equation (eq. 1) one method of solving the time independent equation is to solve the time dependent 
one and let the system relax to its steady state solution. However, we choose a second approach which 
reduces the CPU usage considerably. 
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First, we discretize the volume of the particles, as ui, U 2 , U 3 , .... We arbitrarily choose Vi = i for all i, 
and denote the number of particles of size i by Ui. Doing so, we now observe that the time-dependent 
equation for ni is: 

^-f cr(ui)^/^ni = 0 -/3oniXo + So. (52) 

Note that the 0 term denotes the fact that the particles with the smallest volume cannot be formed through 
the coagulation of smaller particles. For n^+i, the equation is: 

On 1 ^ 

(j{k)^/^nk + (j{k + E mnk+i-i - /3onk+iXo- (53) 

2 = 1 

Since the coagulation coefficient is assumed to be constant, the dynamics of n^+i depends only on 
m ... Tifc and on xo which was found to be ^y2So/Po in §2.2. As a consequence, the equations can be 
solved sequentially from rei. 

Note also that the solution can be accelerated by using FFT to compute the convolution term. In this 
fashion we easily reach ^ 2 , 000,000 in our simulations using just a single CPU. 
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